Band-gap bowing and p-type doping of (Zn, Mg, Be)0 wide-gap semiconductor alloys: 

a first-principles study 



Hongliang Shi 1 !^ and Yifeng Duan 1,2 
1 State Key Laboratory for Superlattices and Micro structures, 
Institute of Semiconductors, Chinese Academy of Sciences, 
P.O. Box 912, Beijing 100083, People's Republic of China 
2 Department of Physics, School of Sciences, China University of 
Mining and Technology, Xuzhou 221008, People's Republic of China 
(Dated: November 28, 2009) 

Using a first-principles band-structure method and a special quasirandom structure (SQS) ap- 
proach, we systematically calculate the band gap bowing parameters and p-type doping properties 
of (Zn, Mg, Be)0 related random ternary and quaternary alloys. We show that the bowing param- 
eters for ZnBeO and MgBeO alloys are large and dependent on composition. This is due to the size 
difference and chemical mismatch between Be and Zn(Mg) atoms. We also demonstrate that adding 
a small amount of Be into MgO reduces the band gap indicating that the bowing parameter is larger 
than the band-gap difference. We select an ideal N atom with lower p atomic energy level as dopant 
to perform p-type doping of ZnBeO and ZnMgBeO alloys. For N doped in ZnBeO alloy, we show 
that the acceptor transition energies become shallower as the number of the nearest neighbor Be 
atoms increases. This is thought to be because of the reduction of p-d repulsion. The No acceptor 
transition energies are deep in the ZnMgBeO quaternary alloy lattice-matched to GaN substrate 
due to the lower valence band maximum. These decrease slightly as there are more nearest neighbor 
Mg atoms surrounding the N dopant. The important natural valence band alignment between ZnO, 
MgO, BeO, ZnBeO, and ZnMgBeO quaternary alloy is also investigated. 

PACS numbers: 71.20.Nr, 71.55.Gs, 61.72.Bb 



I. INTRODUCTION 

Wide band-gap II- VI semiconductors ZnO, MgO, and 
BeO have been extensively investigated because of their 
great potential applications in the spintronic and opto- 
electronic devices [![. Recently ZnBeO alloy has been 
proposed as another wide-band gap oxide, and its band 
gap can be efficiently engineered to vary from the ZnO 
band gap (3.4 eV) to that of BeO (10.6 eV)fl. Further- 
more, ZnO based UV light emitting diodes (LEDs) with 
an active layer region composed of ZnO/ZnBcO quantum 
wells also have been fabricated!?!,!!]. Practically alloying 
is a good approach to perform band-gap engineering to 
extend the available band gap, and the band gap E g of 
the alloy A X B\- X C can be described as 

E g (x) = xE g (AC) + (1 - x)E g {BC) - bx(l - x), (1) 

where b is called the band-gap (optical) bowing param- 
eter, and Eg(AC) and E g (BC) are the band gap of the 
binary constituents AC and BC, respectively. For most 
semiconductor alloys, the bowing parameter b is almost 
independent of composition x^fy. However, the bowing 
parameters b of Gai-^N^As and Zni_ x Te x O alloys are 
strongly composition dependent, indicating that theses 
alloys contain some elements with large differences be- 
tween the size and atomic orbital energies. These alloys 
have attracted a great deal of attention [all]. It is of great 
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interest therefore to investigate how the bowing parame- 
ters b and the band gaps of ZnBeO, MgBeO, and ZnMgO 
to vary since the large differences between Zn (Mg) and 
Be atoms also exist in these alloys. 

It is well known that the p-type doping bottleneck 
poses a challenge for the full utilization of the wide band- 
gap oxides mentioned above. Due to the low valence 
band maximum (VBM) , the transition energies of accep- 
tor defects are always deep and the acceptors are hard 
to ionize at normal operating temperature. Great ef- 
forts have been made to overcome the p-type difficulty, 
for example, through the codopin g ap proach, complex 
defect, and impurity bandfiol fill. Il2||. Li et al. also 
proposed that replacing Zn with isovalent Mg or Be sur- 
rounding with N acceptor defects may act to lower the 
acceptor transition energy through reducing the kinetic 
p-d repulsion In the ZnBeO alloy, it is easy to form 
No+nBez n complexes; thus the transition energy is ex- 
pected to be shallower when the No defect is surrounded 
by more Be neighboring atoms. Therefore it is also in- 
teresting to systematically study the transition energy of 
the No defect in the ZnBeO alloy. 

In this work we systematically investigate ZnBeO, Zn- 
MgO, and MgBeO alloys, in order to understand their 
band structure properties. We hope to show that the 
bowing parameters are large and dependent on compo- 
sition for Mgi-^Be^O and Zni-^Be^O alloys, and much 
larger even in the Be-rich region. However, the bowing 
parameter b is small for the Zni-^Mg^O alloy. We also 
suggest that introducing a small amount of Be into MgO 
will reduce the band gap. We note that these phenomena 
can be explained by the size and atomic orbital energies 



2 



differences between the constituents and isovalent defect 
ievels in these systems. Furthermore, we also show that N 
is a good p-type dopant in the Zno.75Beo.25O alloy which 
may be a good p-type doping material. However, the 
transition energies of No in the Zno.687sMgo.25Beo.o6250 
quaternary alloy lattice-matched to GaN are deep due 
to the lower VBM according to our band offsets calcula- 
tions. 



II. CALCULATIONAL DETAILS AND 
METHODS 

Our first-principles band-structure and total energies 
calculations are performed using the density functional 
theory (DFT) within the local-density approximation (l3l| 
(LDA) as implemented in the Vienna ab initio simula- 
tion package (vasp) code[T3|. The electron and core in- 
teractions are included using the frozen-core projected 
augmented wave approach [15J. The Zn 3d electrons are 
explicitly treated as valence electrons. The whole elec- 
tron wave function is expanded in plane waves up to a 
cutoff energy of 400 eV. All the geometries are optimised 
by minimising the quantum mechanical forces acting on 
the atoms. For the Brillouin zone integration, we use the 
k points that are equivalent to the 4x4x4 Monkhorst- 
Pack special fc-point meshes [l6j in the zinc-blende (ZB) 
Brillouin zone. 

In order to calculate the band structure parameters of 
A X B\_ X C alloys, we adopted the more efficient special 
quasirandom structure (SQS) approach [13, [Hj]- ^ n tms 
approach, we use a smaller unit cell to model the ran- 
dom alloys, where mixed-atom sites are occupied based 
on the physically most relevant structure correction func- 
tions which are closest to the exact values of an ran- 
dom alloy. In our calculation, we use the 2ax2ax2a 64 
atoms SQS at x= 0.0625, 0.25, 0.50, 0.75 and 0.9375 for 
the Mgi-^Be^O, Zni-^Be^O and Zni-^Mg^O ternary 
alloys. 



III. RESULTS 

A. Related binary, ternary and quaternary alloys 

In order to investigate the electronic and structure 
properties of the random alloys, it is essential to cal- 
culate the lattice constants and band gap parameters 
of related ZnO, MgO, and BeO binary alloys. Wei el 
al. concluded that the band structures of the ZB and 
WZ phases for ZnO (GaN) are very similar near the 
band edge at T point through systematically investigat- 
ing the relationships between the band gaps of ZB and 
WZ phase of semiconductors fill] . All our calculations 
are performed for cubic zinc-blende(ZB) alloys and our 
results are also applicable for the ground-state WZ alloys 
similar to other' studies |20j|. 



TABLE I: Calculated lattice constant a, bulk modulus B, and 
band gap E g at F for ZB binary alloys ZnO, MgO, and BeO. 
The experimental values of B, and E g suggested in Ref. 7 are 
also listed in parentheses. 





ZnO 


MgO 


BeO 


a(A) 


4.506(4.47°) 


4.517(4.524") 


3.766(3.768 c ) 


B(Mbar) 


1.562(1.837) 


1.554(1.603) 


2.354(2.244) 








2.29 c 


E g (eV) 


0.70(3.4) 


3.731(7.67) 


7.852(10.585) 



"experimental value in Ref. 20 
^theoretical value in Ref. 21 theoretical value in Ref. 22 



Table I shows our calculated equilibrium lattice con- 
stant a, bulk modulus B, and band gaps at T point for 
zinc-blende structure ZnO, MgO, BeO. We also list the 
corresponding experimental values of bulk modulus B 
and band gaps (in brackets) for their ground-state struc- 
ture. Our calculated lattice parameters for aznO, flMgOi 
and aecO are 4.506, 4.517, and 3.766 A, respectively, 
which are in good agreement with the experimental and 
theoretical values 4.47, 4.524, and 3.768 A. 

Conventionally, the lattice constant of Zni^Be^O al- 
loy is described by the Vegard law(23|, that is 

a(x) = (1 - x)a Zn o + xa BcQ . (2) 

Note that the lattice mismatch between ZnO (MgO) and 
BeO is 19.6% (19.9%). Consequently, the ZnBeO and 
MgBeO alloys are highly strained and the lattice con- 
stants may deviate from the Vegard law. First, we calcu- 
late the lattice constant by the Vegard law, then fit the 
Murnaghan equation of state [24[ at various lattice con- 
stants around it. Figure 1 shows the lattice constant as 
a function of concentration x obtained from the Vegard 
law (black dots) and fitted by the Murnaghan equation 
(red dots). We note that (i) For the ZnBeO and Mg- 
BeO alloys, due to the large size difference between the 
Zn (Mg) and Be atom, the lattice constants calculated 
by LDA (by fitting the Murnaghan equation) are slightly 
larger than those obtained from the Vegard law. This 
is because the Zn-0 (Mg-O) bond is hard to compress. 
Our calculated results also show that nearly all the Be-0 
bonds are longer than those of Be-0 in the host BeO. 
The biggest deviation occurring at composition x=0.5 is 
2.5% for ZnBeO alloys, (ii) the Vegard law is a good 
model for ZnMgO alloys due to the lattice match and 
the similar size of Zn and Mg ions. In actuality, we use 
the lattice constants obtained by fitting the Murnaghan 
equation in all our calculations for ternary alloys. 



B. The bowing parameters b of ZnMgO, ZnBeO, 
and MgBeO alloys 

We calculate the LDA band gaps to be 0.70, 3.731, 
7.852 eV at T point for zinc-blende ZnO, MgO and 
BeO, respectively. Due to the LDA underestimation of 
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FIG. 1: Comparison between Vegard's law and LDA calcu- 
lated (by fitting the Murnaghan equation of state) lattice con- 
stants of Mgi-^Be^O, Zni-zBe^O and Zni-zMg^O alloys. 
The picture on the left and right corresponds to two endpoint 
alloys (a) MgO (x=0) and BeO(z=l), (b) ZnO (x=0) and 
BeO(a;=l), and (c) ZnO (x=0) and MgO(:r=l), respectively. 
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FIG. 2: The predicted band gaps of Mgi-^Be^O, Zni-^Be^O 
and Zni-zMgzO ternary alloys calculated at various compo- 
sitions. 



TABLE II: Calculated band-gap bowing coefficient b (eV) for 
ternary alloys at different compositions. 

composition x 

Ternary Alloys 0.0625 025 O50 075 0.9375 

Zm_ I Mg I 0.8794 0.9797 1.4172 2.2827 4.3430 
Zm-^Be^O 4.8215 5.8189 8.2116 12.8925 23.9648 

Mgi^Be^O 4.4951 5.2883 5.4772 7.2861 18.0994 



band gap, our values are smaller than the experimental 
valuesQ of 3.4, 7.67, 10.585 eV for their ground-state 
phase, respectively. However, the error of bowing pa- 
rameter b is small in Eq. 1 because the systematic band 
gap error cancels out in deriving the bowing parameters 
through comparing chemically identical systems in two 
different forms (AX and BX1 [2al . This has been con- 
firmed by numerous studies [{j l2aj. 

Table II presents the calculated band-gap bowing co- 
efficients for Mgi-^Be^O, Zni_ x Be x O and Zni_ x Mg x O 
ternary alloys as a function of the composition x. our 
results show some chemical trends. 

(i) The bowing parameters b are large for ZnBeO and 
MgBeO alloys. This is thought to be because of the 
atomic size difference and the large chemical mismatch 
between Zn (Mg) and Be. However, compared with Zn- 
BeO and MgBeO the bowing parameters of ZnMgO al- 
loys are smaller as Zn and Mg are of similar size. The 
atomic size of Zn is a little larger than that of Mg, thus 
the bowing parameters of ZnBeO are also a little larger 
that those of MgBeO. From Be to Mg to Zn, the atomic 
size difference increases as does the atomic p-orbital en- 



ergy difference. 

(ii) The bowing parameters are composition and large 
at Be-rich region for ZnBeO and MgBeO alloys. Taking 
ZnBeO as an example, the bowing parameter b is 4.8215 
eV in the Zn-rich region, whereas in the Be-rich region 
the b is 23.9648 eV. This is because the size and chemi- 
cal differences between Zn and Be are large. Incorpora- 
tion of a small amount of Zn into BeO at substitutional 
sites creates a deep isovalent defect level inside the band 
gap. Our direct calculations show that the position of 
the isovalent defect level is at 0.66 eV below the conduc- 
tion band minimum (CBM). Furthermore the VBM is 
pushed up due to the p-d repulsion introduced by Zn 3d 
orbital. This leads to a large band-gap reduction and a 
large bowing parameter in the Be-rich region. However, 
upon adding a small amount of Be into ZnO, the VBM 
will shift downwards by a small amount due to the reduc- 
tion in p-d repulsion. In this case the band gap increases 
a little and the bowing parameter is small. 

(iii) Using the bowing parameters and experimental 
band gaps, we calculate the band gaps of ZnBeO, Mg- 
BeO, and MgBeO alloys at different concentrations. The 
results are summarized in Figure 2. We notice that al- 
though the band gap of BeO is larger than that of MgO, 
incorporation of a small amount of Be into the MgO ma- 
trix will lead to a reduction in the band gap, indicating 
that the bowing parameter is larger than the band-gap 
difference. For Mgo.75Beo.25O, the bowing parameter b 
(5.2883 eV) is larger that the band gap difference (2.915 
eV). 
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C. The formation energies and transition energy 
levels of N in ZnBeO and ZnMgBeO quaternary 
alloys 

P-type doping of ZnO is difficult due to the low VBM 
resulting from the largely electronegative characteristic of 
oxygen. In order to lower the acceptor ionization energy 
extensive studies have been performed. Park et al.\lty 
showed that N is the best p-type dopant source among 
the group-I and V elements. Li et aL[ll[ also showed 
that replacing Zn with Mg or Be surrounding the No 
substitutional defect can decrease the acceptor transi- 
tion level. In the ZnBeO ternary alloy, it is easy to form 
No+nBezn complexes; thus we expect the acceptor tran- 
sition level becomes shallower when the No is surrounded 
by more Be atoms. In the following study, we choose the 
Zno.75Bco.25O alloy approximately as the experimental 
investigations do[5[, and we also choose the N as dopant 
due to its low valence p-orbital energy. 

The ionization energy of an accepter (q<0) with re- 
spect to the VBM in the impurity limit is calculated by 
the following procedure, further described in Refs. 26 
and 27, 

e(0/q) = {E(a, q) - [E(a, 0) - qe k D (0)]}/(-q) 

+ [^>(0)-e vBM (host)], (3) 

where E(a,q) and E(a,0) are the total energies of the 
supercell at charge state q and neutral, respectively, for 
defect a, e h D (Q) and e^(0) are the defect levels at the spe- 
cial k points (averaged) and at the V point, respectively, 
and ey BM (host) i s the VBM energy of the host supercell 
at the r point and is aligned using core electron levels 
away from the defect. The first term on the right-hand 
side of Eq. (5) determines the U energy parameter (in- 
cluding both the Coulomb contribution and atomic re- 
laxation contribution) of the charged defects calculated 
at the special k points, which is the extra cost of en- 
ergy after moving (-q) charge from the VBM of the host 
to the neutral defect level. The second term gives the 
single-electron defect level at the T point. For charged 
defects, a uniform background charge is added to retain 
the global charge neutrality of the periodic unit-cell. 

The defect formation energy of a neutral defect is de- 
fined as 

AH(a, 0) = E(a, 0) - S(host) + n Mo + n Q /i Q , (4) 

where E(host) is the total energy of the host supercell 
without defect a; /zo and \i a are the chemical potentials 
of constituents O and a (N) relative to the element gas 
energy. The no (no>0) and n a (n Q <0) are the numbers 
of O and extrinsic defects a (N). We calculate the chem- 
ical potentials for O(^o) and N(^n), relative to gaseous 
O2 and N 2 . 

Table III shows our calculated No acceptor formation 
energies and e(0/-) transition energy levels for the ZnO 
and Zno.75Beo.25O alloy. For N in Zno.75Beo.25O alloy, 
we calculate the Nq at different sites. The nNN in Table 



TABLE III: Calculated formation energy AH/(«,0) (eV) and 
transition energy level e(0/-) (eV) for Zno.7sBeo.2sO:N and 
ZnO:N are listed. 





AH/(a,0) 


e(0/-) 


ONN 


4.793 


0.533 


INN 


4.818 


0.347 


2NN 


5.134 


0.233 


ZnO 


4.661 


0.362 



III denotes that there are n Be atoms in the tetrahc- 
dral nearest neighbor (NN) sites centered around the N 
atom. Our calculated results show the following typical 
behaviours. 

(i) The formation energy of neutral No defect in the 
Zno.75Beo.25O alloy increases when there are larger num- 
bers of Be atoms in the defect neighboring sites. This 
is because the No substitution introduces a larger strain 
into the Be-rich region than into the Zn-rich region. The 
formation energy of No in the ZnO is smaller than that 
in the alloy. This is because the volume of ZnO is about 
0.9% greater than that of the Zno.75Beo.25O alloy, in- 
dicating that an No defect in pure ZnO introduces a 
smaller strain. 

(ii) For No in ZnO, our calculated No e(0/-) acceptor 
level is at 0.362 eV above the VBM, which ag rees well 
with previous first-principles calculations (TollTlj|. For the 
case of the Zno.75Beo.25O alloy, we note that the position 
of No e(0/-) acceptor level is dependent on the site of the 
defect. It decreases as the number of N nearest neigh- 
bouring atoms increases, from 0.533 eV when N has no 
nearest neighboring (ONN) Be atoms to 0.233 eV when 
there are 2 NN Be atoms. This is thought to be because 
the acceptor level is lowered by the reduction in the p-d 
repulsion. A more detailed explanation is presented. For 
ZnO, the VBM state consists mostly of po, PZn, and dzn 
orbitals. For an No acceptor in ZnO, the defect level 
mainly consists of N impurity valence p orbitals, and the 
acceptor level is pushed up by the p-d (pN-dzn) repul- 
sion introduced by the same symmetry. However, for Be, 
it has no occupied 2p and 3d orbitals. If Zn is substi- 
tuted by Be, the p-d (pN-^Zn) repulsion will decreases 
and the No acceptor transition level becomes shallow. 
We also note that the p-p repulsion can push the accep- 
tor level further down, although it is weaker compared 
with p-d repulsion because the p-p (p-d) repulsion is in- 
versely proportional to the energy difference between the 
Pn and ps c (pn and rfzn) valence orbitals[28]. Similarly, 
in the Zno.75Beo.25O alloy, from ONN to 2NN, the p N - 
dzn repulsion decreases and the p-p (p^-PBe) repulsion 
increases; thus, the No acceptor level becomes increas- 
ingly shallower. In ONN, the No acceptor level in the 
Zno.75Beo.25O alloy is deeper than that in ZnO due to 
the lower VBM (see below). 

(iii) According to Eq. 4, we see that the chemi- 
cal potentials have large effect on the defect formation 
energy. Practically, the chemical potentials can vary 
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in a range because of some thermodynamic limitsf^. 
The high defect formation energy leads to a low dopant 
solubility, which may be overcome by non-equilibrium 
growth methods[29]. If the formation energy of No in 
Zno.75Beo.25O alloy is reduced so that it is smaller than 
that in ZnO by controlling the chemical potentials of N 
and O during the growth process, there are more N atoms 
residing in the INN and 2NN sites with low acceptor 
transition levels. Based on this, we predict that ZnBeO 
alloy may be a better p-type doping material than ZnO. 

We have also investigated the electronic structures and 
p-type doping properties of the ZnMgBeO quaternary al- 
loy lattice-matched to GaN substrate. For quaternary 
alloys, the band gap and lattice constant can be indi- 
vidually tuned. For Zni-^-yMgyBe^O alloy, the lattice 
constant can be described by 



a(x, y) = (1 - x - y)a Zn0 + ya Mg o + zaBcO ■ 



(5) 



For the Zn22MggBe2032 alloy, we find the Vegard law is 
well obeyed. We also note that the lattice constants of 
Zn22Mg 8 Be 2 032 quaternary alloy and zinc-blende GaN 
are 4.461 and 4.462 A[3(J, respectively, thus the lattices 
match well. 

The band gap of the Zni-^-yMgyBe^O quaternary al- 
loys is a function of compositions x and y and can be 
described by a second order expansion in x and y[3lj. 



E g (x,y) = (l-x-y)E g (ZnO) + yE g {MgO) 
+xE g (BeO) - &z„M g oy(l - y) 
-fenBcO^(l - x) - b xy xy. 



(6) 



In order to get b xy approximately, we use the SQS at 
(z,y)=(0.0,G.5), (0.5,0.0), and (0.5,0.5). The parameter 
b xy is obtained using b xy =b Zn B C o-bznMgO-bMgBoO under 
a quadratic approximation. Our calculated b xy is 1.317 
eV. Using the experimental band gaps and Eq. 6, we 
predict that the band gap of Zn22MggBe2032 alloy is 
4.149 eV. 

Table IV shows calculated No acceptor forma- 
tion energies and e(0/-) transition energy levels for 
Zn 2 2Mg 8 Be 2 032 alloy. The nNN in Table IV denotes 
that there are n Mg atoms in the tetrahedral near- 
est neighbor (NN) sites centered around the N atom. 
The ONNlBe denotes that N atom has one Be atom 
and without Mg atoms in the NN sites. From ONN to 
3NN, the e(0/-) decreases slightly less than that of the 
Zno.75Beo.25O alloy. This is because the pm? orbital en- 
ergy level is higher than that of PBe[28|], thus the p-p 
repulsion pushing down the acceptor level is weaker in 
the quaternary alloy. The e(0/-) is smaller in ONN 1 Be 
than that in ONN, which is also consistent with our find- 
ings for No in the Zno.75Beo.25O alloy. We note that the 
No e(0/-) transition energy level is deeper than that in 
ZnO because of the large downward shift of the VBM 
according to our band alignment calculation (see below) . 
The CMB of Zn 2 2Mg 8 Be 2 032 alloy is as low as that of 
ZnO, indicating that it is also can be n-type doped easily. 



TABLE IV: Calculated formation energy AH/(a,0) (eV) and 
transition energy level e(0/-) (eV) for ZnMgBeO:N are listed. 





AHj(q,0) 


e(0/-) 


ONN 


4.773 


0.546 


INN 


5.176 


0.489 


2NN 


5.611 


0.458 


3NN 


6.033 


0.446 


ONNlBe 


4.978 


0.477 



IV. THE NATURAL VALENCE BAND 
ALIGNMENT 

Band offset is an important parameter in the het- 
erostructures. We shall describe how we investigated the 
natural band alignment by adopting the approach pro- 
posed by Wei [33]. This method is also widely used in 
photoemission core level spectroscopy [33J . According to 
Wei's method [32|. the band offset is given by 



AE V (AX/BY) = A£ By , - AE, 



? AX 
J v,C 



Here. 



E 



AX 



E, 



AX 



AE* X J BY . (7) 



(8) 



(and similarly for AE^, ) are the energy separations be- 
tween core level (C) and valence band maximum energy 
for pure AX (and similarly for pure BY), while 



A K X J BY = E c Y - E c X 



(9) 



is the difference in core level binding energy between AX 
and BY at the AX/BY superlattice. We calculate the 
AE^q and AE^y, at the equilibrium zinc-blende lattice 
constants for AX and BY, respectively. We use the aver- 
age lattice constant a (between AX and BY) to calculate 
the AE^ X J, BY in the AX/BY superlattice. A more de- 
tailed description of the procedure can be found in Ref. 
32. 

The natural conduction band minimum (CBM) align- 
ment is calculated by 

AE C (AX / BY) = AE V (AX / BY) + AE g (AX /BY), (10) 

where AE g (AX/ BY) is the experimental band gap dif- 
ference between the AX and BY alloys. 

Our calculated natural band alignment is shown in 
Fig. 3. For the VBM alignment, the VBM of ZnO is 
higher than that of MgO and BeO. This is because the 
p-d {po — G?Zn) repulsion pushes up the VBM, whereas 
Mg and Be have no active d electrons. According to the 
common-anion rule[32], the valence band offset is small 
between common-anion systems. However, the VBM of 
BeO is higher than that of MgO. This is because BeO 
shows strong covalency, resulting in the kinetic-energy- 
induced valence band broadening 34] . If we assume that 



6 




MgO BeO Zn 22 Mg 8 Be 2 Zn 24 Be a O ZnO 



FIG. 3: The calculated natural band alignment of MgO, BeO, 
Zn22Mg8Be2 032, Zn24Bes032 and ZnO alloys. 

the VBM varies linearly as a function of the alloy com- 
position, Zn24Be§032 would have a VBM approximately 
0.32 eV lower than that of ZnO. However, our direct cal- 
culation shows that VBM of Zn24Be§032 (Zno.75Beo.25O) 
alloy is only about 0.11 eV lower than that of ZnO due 
to the large p-d repulsion. In order to investigate the 
natural conduction band alignment, we use the experi- 
mental band gaps demonstrated in Fig. 3. We see that 
the Zno.75Beo.25O alloy acts as a good barrier material for 
ZnO-based UV LEDs. The AE C and AE V have values 
of approximately 0.60 and 0.11 eV, respectively, which 
are large enough to strongly confine the electrons and 
holes well|35|. The quaternary alloy Zn22MgsBe2032 has 
a CBM as low as that of ZnO, indicating that it can be 
n-type doped easily. By contrast p-type doping is more 
difficult due to the low VBM which is also confirmed by 
our study of p-type doping above. Based on these chemi- 



cal behaviours, we can tailor the band offset to practically 
perform the band-engineering as desired. We hope that 
our results will be helpful for the applications of opto- 
electronic devices. 



V. SUMMARY 

We have investigated the band gap bowing parame- 
ters and p-type doping properties of (Zn, Mg, Be)0 re- 
lated random ternary and quaternary alloys using first- 
principles band-structure methods. The bowing param- 
eters are large and dependent on composition for ZnBcO 
and MgBeO alloys due to their large size and chemical 
mismatch. We also show that the band gap of MgO 
will reduce on adding a small amount of Be into MgO. 
The acceptor transition energy level becomes shallower 
when N has more Be nearest neighbors, thus ZnBeO 
(Zno.75Beo.25O) may be a good p-type material. The 
electronic and p-type doping properties of ZnMgBeO 
quaternary alloy lattice-matched to GaN substrate are 
also studied. According to our natural band alignment 
calculations, n-type doping of Zn22MggBe2032 should be 
easy due to the low CBM. We hope our conclusions have 
applications to optoelectronic devices. 
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